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Summary. With polarization of the cosmic microwave background (CMB) now 
detected, and confirmed by several independent experiments, the next goal is to 
characterise accurately its statistical properties. In these lecture notes we review 
the physical motivation for pursuing CMB polarization, and the basic statistical 
properties of the polarization fields. We then discuss some of the key aspects of the 
analysis of CMB polarization data, focusing on the additional complications that 
arise compared to temperature data due to the tensor character of the polarization 
field. 



1 Introduction 

The observed statistical properties of the temperature anisotropies in the cos- 
mic microwave background (CMB) have played a key role in shaping our un- 
derstanding of the large-scale properties of the Universe. The angular power 
spectrum of the anisotropies, quantifying the contribution to the variance 
of the anisotropies as a function of angular scale, has now been measured 
exquisitely over two decades of scale by the Wilkinson Microwave Anisotropy 
Probe (WMAP) satellite [Q. A number of ground and balloon-borne experi- 
ments extend the power spectrum measurements by a further decade in scale 
but with somewhat less precision than has been attained on the larger scales 
(see for a recent review) . Current CMB data is remarkably consistent with 
predictions based on the passive evolution of density perturbations, with adi- 
abatic initial conditions drawn from a Gaussian distribution with an almost- 
scale-invariant power spectrum, in a spatially-flat universe. Such primordial 
fluctuations, and the observed flat geometry of space, are natural outcomes 
of the simplest models of inflation in the early universe, e.g. 

In addition to having angular variations in the total intensity, with r.m.s. 
~ 120 /xK, the CMB is also partially linearly polarized with r.m.s. ~ 6 //K. 
Polarization is generated by Thomson scattering of CMB photons around 
the time of recombination, once scattering becomes sufficiently rare to allow 
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the development of anisotropics in the CMB intensity. Polarization is further 
generated at late times by re-scattering once the Universe reionizes. Despite 
being predicted shortly after the discovery of the CMB 4 , linear polarization 
wasn't detected until 2002 by the Degree Angular Scale Interferometer (DASI) 
team |Kj . It has now been measured by two further groups |HJ [7] , and also de- 
tected via its correlation with the temperature anisotropies in the first-year 
WMAP data [BJ- The direct measurements are currently much less precise 
than those for the temperature anisotropies, but they are fully consistent 
with predictions in models favoured by the temperature data. Future, precise 
measurements of CMB polarization will be very valuable, allowing accurate 
determination of several cosmological parameters that are almost degener- 
ate in their effect on the temperature anisotropies [H]. The biggest winners 
will be the reionization history |10| . with the temperature-polarization cor- 
relation observed by WMAP already constraining the optical depth back to 
reionization to be ~ 0.15 [5], and constraints on any gravitational- wave back- 
ground ^2 1X3 and sub-dominant isocurvature modes |13| . 

As the accuracy and size of CMB datasets have grown, so have the de- 
mands on the techniques used to process the data. A standard set of com- 
pression and cleaning steps have emerg ed that ultimately distill O(10 10 ) time- 
ordered observations to several hundred power spectrum measurements, and 
only around 10 cosmological parameters and their errors (see e.g. ^3] for a 
review.) Although the detail of some of these steps are instrument-specific, 
most analyses can be regarded as (approximate) variants on a small number 
of generic algorithms. Three of the most important steps are reviewed in de- 
tail elsewhere in this volume: for map making and power spectrum estimation 
see the contribution by Borrill; for (astrophysical) component separation and 
foreground removal see that by Delabrouille. Most of the steps in the analysis 
of CMB polarization data are straightforward generalisations of those for the 
temperature, but some new issues do arise mainly due to the tensorial nature 
of the polarization field. Furthermore, the small amplitude of the polariza- 
tion signal demands that even more careful attention be paid to potential 
systematic effects in the analysis of polarization data than for the temper- 
ature anisotropies. Inevitably, this makes polarization data processing more 
instrument-specific and was the reason that the WMAP team did not release 
polarization maps with their first-year data. 

The purpose of these lecture notes is to describe, in outline, the scien- 
tific motivation for pursuing CMB polarization, and discuss some of the key 
steps in the analysis of polarization data. We shall focus on those parts of the 
analysis that are complicated by the different geometric character of linear 
polarization compared to temperature anisotropies. To keep the discussion 
general, we will not have much to say on detailed instrument-specific parts 
of the analysis, although these are clearly crucial to the success of any given 
observation. We start in Sect.[5]with a discussion of the statistics of the CMB 
polarization fields and their representation in terms of orthonormal basis func- 
tions on the sphere. Section previews the basic physics of CMB polarization, 
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and the current detections and best upper limits on CMB polarization are 
discussed in Sect. 0] What we have already learnt from current polarization 
measurements, and what we may learn in the future, is described in Sect. |SJ 
The analysis of polarization data is then treated in Sect. El where discussions 
of map-making, astrophysical component separation, E- and -B-mode sepa- 
ration, power spectrum estimation and non-Gaussian lensing issues can be 
found. 

2 Statistics of CMB Polarization 

The polarization on the sky along some line of sight ft is conveniently described 
in terms of Stokes specific brightness parameters I(n), Q[n), U(h) and V(n). 
For the CMB the total intensity I(n) has a Planck spectrum along any line 
of sight, but with angular temperature variations AT(n) at the 10 -5 level on 
a 2.725-K background. In linear theory the anisotropic contribution to I(n) 
is thus proportional to AT(n) and has a spectrum given by the derivative 
of the Planck function (with respect to temperature) . The parameters Q and 
U describe linear polarization, and inherit the spectrum of the anisotropic 
part of I(n) since CMB polarization is produced by frequency-independent 
Thomson scattering. The parameter V describes circular polarization and is 
expected to vanish for the CMB. 

For a line of sight n, the radiation is propagating along —n. If we introduce 
a pair of orthogonal directions in the surface of the sphere, and call these 
x and y, such that x, y and —h form a right-handed basis, then the linear 
Stokes parameter Q can be defined operationally as the difference between the 
intensities transmitted by perfect polarizers oriented along x and y. Similarly, 
U is the difference in intensities if the polarizers are rotated by 45°. We shall 
adopt the convention that in a spherical-polar coordinate system, the local 
x-axis is along the direction of decreasing polar angle 9 (i.e. north) and the 
y-axis is along the direction of increasing azimuthal angle cf> (i.e. east). This is 
in accordance with the IAU recommendations, but differs from e.g. and 
the latest version (1.2) of the HEALPix package 1 . If the local x- and y-axes 
are rotated through an angle ip in a right-handed sense about the propagation 
direction, Q and U transform as 



These show that Q and U properly form the components (in an orthonormal 
basis) of a rank-2, symmetric trace-free tensor V (n) that is transverse to 
the line of sight: 



Q -> Q cos 2ip + U sin 2ip , 
U -> U cos 2ip - <3sin2V> . 



(1) 




(2) 



i 



http://www.eso.org/science/healpix/ 
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Fig. 1. Polarization patterns for a pure-electric mode (left) and pure- magnetic mode 
(right) on a small patch of the sky for potentials that are locally Fourier modes. 
The shading denotes the amplitude of the potential. For the electric pattern the 
polarization is aligned with or perpendicular to the Fourier wavevector depending 
on the sign of the potential; for the magnetic pattern the polarization is at 45°. 

This linear polarization tensor is proportional to the correlation tensor (at 
zero lag) of the electric field components for the radiation propagating along 
— n. 

Any two-dimensional tensor with the symmetries of V ab can be derived 
from two scalar fields, denoted Pe and Pb, via [111 112) : 

Vab = V <a V b) P B + € c (a V b) V c P B , (3) 

where angle brackets denote the symmetric, trace-free part of the enclosed 
indices, V a is the covariant derivative on the sphere, and e a b is the alternating 
tensor. The part derived from Pe is the electric (or gradient) part and that 
from Pb the magnetic (or curl) part of V a b- The divergence V a V a b is a pure 
gradient if Pb = 0, and a curl if Pe — 0. The decomposition Q is analogous 
to writing a vector field as a gradient and a divergence free vector, since in two 
dimensions the latter is always of the form e a b ^7bX- To gain some intuition for 
E modes and B modes, consider the case where Pe and P B behave locally like 
a plane wave across a small patch of the sky (that can accurately be treated as 
fiat). The electric and magnetic contributions to the polarization are depicted 
in Fig. ^ Quite generally, for a given potential P, the magnetic polarization 
V^ b = e c ^ a V f) )V c P generated from P is related to the electric polarization 
= V (a V b) P by Vg = e\ a V^ )c which is a right-handed rotation of Vg by 
45° about the line of sight. 

The fields Pe and P B are simple scalars that can be expanded in spherical 
harmonics with Z > 2 as 




Bi m Yi m {h) . 



(4) 

(The normalisation is conventional.) The polarization tensor can then be writ- 
ten as 

Vabin) = ±= Ei m Y { f m)ab + B lm Y ( f m)ab , (5) 

Irn 

where the trace-free, symmetric tensors 
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E 2(1-2)1 

r (hn)ab = y (/_|_2)! V ( aVb ) rim 



Y (f,n)ab = y y + 2 ^,' e C (aV b )V e Yi m , (6) 

form orthonormal tensor bases for electric and magnetic polarization respec- 
tively. The electric and magnetic multipolcs, E\ m and Bi m , can be extracted 
directly from the polarization tensor by contraction with the appropriate ten- 
sor harmonic followed by integration over the sphere. For example, 

E lm =V2jr ab Y { fZ )ab dn. (7) 

Methods of recovering E and B modes when observations cover only a limited 
part of the sky are discussed later in Sect. El Reality of Pe and Pb demands 
that E* m = (-l) m Ei _ m , and similarly for Bi m . 

We end this subsection by noting an alternative formalism for expressing 
symmetric, trace- free tensors on the sphere: the spin- weighted formalism |15j . 
first introduced to CMB physics in The complex polarization P = Q + iU 
is scaled by a phase factor e~ 21 ^ under a right-handed rotation of the local 
x- and y-directions about the propagation direction, and so is defined to have 
spin -2. The appropriate basis functions for expanding functions of definite 
spin on the sphere are the spin- weighted spherical harmonics, ±2Yi m (fi), which 
are related to the components of the tensor harmonics, Y^ m ^ ab and Yff m y ab , 
in a null diad: 

1 

i 

Here the null (complex) vector m = (8 + i(j))/^/2 satisfies m ■ m = and 
m ■ m* = 1. With Q and U expressed on the north-east basis of a spherical- 
polar coordinate system, the multipole expansion of the complex polarization 
is 

(Q ± iU)(n) = J2( E im T iB lm ) T2 Yi m (n) . (9) 

lm 

The electric and magnetic multipoles can be extracted using the orthonormal- 
ity of the spin-weighted harmonics, e.g. 

Ei m T iBi m = J(Q ± iU) T2 Y t * m dn . (10) 
2.1 Polarization Power Spectra 



Y (f m )ab = -^{-2Yi m m a m b + 2 Yi m m* a m* b ) 
Y (? m )ab = T-7^(-2Yi m m a m b - 2Y hn m* a m* b ) . (8) 



The decomposition of the polarization field into electric and magnetic parts 
is invariant under rotations, and the electric and magnetic multipoles at fixed 
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I transform irreducibly under rotations, i.e. Ei m — > D l mm ,Ei m i where D l mm , 
are the Wigner functions (see e.g. JHI)- Under the operation of parity, (Q ± 
iU)(h) — > (Q =F ill)(—n) so that E lm — > (— l) l Ei m (electric parity) while 
Bi m — » (— l) l+1 Bi m (magnetic parity). These transformations ensure that 
the potential Pe is a scalar under parity, Pe(n) — > Pe(— n), but Pg is a 
pseudo-scalar, Pe(n) — > — Pb(— w). The temperature anisotropies AT(n) are 
a scalar function and so can be expanded in spherical harmonics in the usual 
way: AT(n) = $Zj m Ti m Yi m . The multipoles T; m have electric parity. 

In the absence of parity-violating interactions, the homogeneous and 
isotropic background universe on which cosmological perturbations are as- 
sumed to propagate can support an ensemble of fluctuations that is invariant 
under parity and rotations. What this means is that in the ensemble any 
fluctuation is just as likely as its parity-reversed, or rotated counterpart. The 
implication of this for the statistics of the CMB is that any expectation value 
must be preserved if we rotate or invert the fields involved. For the two-point 
statistics, this limits the non-zero correlations between the observable multi- 
poles to 

(TlmT^^) = 5u'5 mm iCf , (El m E*, m ,) = Su'S mm 'C E , 
(Bi m Bi, m ,) — Sll'Smm'Cj 3 , {Ti m Ei, m ,) — 5ll'8mm'Cf E . (If) 

For Gaussian fluctuations, such as those predicted by the simplest models of 
inflation, these power spectra Cf, Cf , Cf and Cj E contain all of the statis- 
tical information in the CMB. There have been several claims of violations of 
Gaussianity/rotational invariance in the one- year WMAP temperature data, 
e.g. |171 1181 H^l 1201 12f| . but the issue of whether these are of primordial ori- 
gin, or rather artifacts of the instrument or local Universe, is unresolved |22j . 
Example power spectra of the CMB observables for adiabatic density pertur- 
bations in a flat, ylCDM model are shown in Fig. The physics that gives 
rise to these spectra is reviewed briefly in Sect. [3] 

It is also interesting to consider the two-point correlations between the 
Stokes parameters (and also the temperature anisotropies) in real space \2',\\ . 
For polarization, the correlation will be manifestly rotationally-invariant only 
if we work with Stokes parameters defined on a basis intrinsic to the two 
points, fi\ and fi2, under consideration. If we define the local x-direction by 
the geodesic connecting the two points, and denote the polarization on these 
bases by an overbar, the correlations are [T^l |2~H l2"5j 

(Q(hMn 2 )) = i^^[Cf(4 2 + 4- 2 )(/9) + Q B (4 2 -4- 2 )(/3)] , 

(U^pinz)) = \J2 ^-i c i E ( d 22 4-a)G3) + Q B (4 2 + 4 -2)^)] . 
/ 

(T(rH)Q(n 2 )) = £ ^± C f E ^M , 
1 
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(Tin^Tin,)) = £ -^±Cfd l M , (12) 

i 

where (3 is the angle between the two points (0 < (3 < ir), and d l mn ((3) are 
the reduced Wigner functions. The correlations (QU) and (TU) vanish due 
to parity and rotational invariance. Examples of the non-zero polarization 
correlation functions are plotted in Fig. |3J along with a simulated Gaussian 
realisation of the temperature anisotropics on a flat patch of the sky with the 
correlated part of the polarization overlaid. As first noted in |23) . the nega- 
tive tail to the temperature-polarization correlation function on large scales 
is generic, arising from the infall of the photon-baryon plasma into poten- 
tial wells around recombination, and gives a tangential pattern of correlated 
polarization around large-scale temperature hot spots. On smaller scales the 
sign of the cross-correlation oscillates. For adiabatic perturbations on scales 
smaller than ~ 0.3°, the cross-correlation is positive and this gives rise to the 
tangential polarization pattern around small-scale cold spots that can be seen 
in Fig. 

Cosmic Variance 

Since for Gaussian CMB fluctuations all cosmological information is encoded 
in the various power spectra, the precision with which the CMB can constrain 
cosmology is determined, in part, by the accuracy to which these spectra can 
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Fig. 3. Polarization correlation functions for adiabatic density perturbations (left). 
The right-hand plot overlays a simulation of the correlated part of the polarization 
on the CMB temperature anisotropics. 



be determined. In the absence of instrumental noise, and with full sky cov- 
erage, cosmic variance provides a fundamental limit to this accuracy. Cosmic 
variance arises from having a single realisation of a non-ergodic process (the 
sphere is compact) from which to estimate the power spectra. In this case, 
the obvious estimators of the power spectra, e.g. 



C E - — 



^yEi^i 2 ' ( 13 ) 



are optimal in that they have the minimum variance of all unbiased estimators. 
At each I we have only 21 + 1 independent real numbers from which to estimate 
a variance, so the cosmic variance in our estimator Cf for Gaussian fields is 

var((7f ) = ^{Cf f . (14) 

Equivalent results apply to the two other auto-spectra Cf and Cf . For the 
cross-spectrum, the cosmic variance is only slightly more complicated |12II26| : 

var( ^ T£) = 27TT [{cIEf + cJcf ] ■ (15) 

The second term arises from chance correlations in the single realisation that 
we have available, and would fundamentally limit the accuracy of searches 
for parity violations via e.g. the cross-correlation between the temperature 
anisotropies and magnetic polarization. Due to the presence of correlations 
between the temperature and the electric polarization, estimates of Cj and 
Cf have non-vanishing covariance, and each is also correlated with the esti- 
mator for Cf E : 
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cov(Cf , Cf ) = ^(Cr) 2 , cov ( cf^ E , Cj E ) = ^A_cf /B CP . 

(16) 

Of course, in the realistic case of non-zero instrument noise and limited sky 
coverage the precision of power spectrum estimates falls short of the cosmic 
variance limit. The effect of only observing a fraction f^y of the sky is not only 
to increase the variance [27], but also to correlate estimates across Al given 
roughly by the inverse linear size of the retained portion of the sky |281 . 
The increase in variance can be crudely accounted for by reducing the number 
of independent modes to (2Z + 1)/ S k y . For polarization, this mode-counting ar- 
gument ignores the loss of ambiguous modes that cannot be properly classified 
as electric or magnetic with partial sky coverage (see Sect. EJ- This effect is 
minimised for surveys covering large connected regions of the sky. 



3 Physics of CMB Polarization 

When the mean free path was small compared to the spatial scale (i.e. wave- 
length for a plane wave) of the perturbations, Thomson scattering kept the 
CMB radiation isotropic in the rest frame of the electron-baryon plasma. 
As protons and electrons started to recombine to form neutral hydrogen, the 
mean free path increased and anisotropies started to develop. Subsequent 
Thomson scattering of the quadrupole component of the anisotropy generated 
linear polarization; this was then preserved for the free-streaming CMB until 
the Universe was reionized by the ionizing radiation from the first non-linear 
structures. The photon visibility function - the probability that a photon last 
scattered around recombination as a function of cosmic time peaks around 
370 kyr after the big bang and has a width ~ 115 kyr |30j. Since polarization 
is only generated by scattering it provides a very clean probe of conditions 
around the time of recombination and reionization. 

Consider unpolarized radiation with a temperature distribution 0(e), 
where e is the radiation propagation direction. If the radiation field pos- 
sesses a quadrupole anisotropy ©2(e), Thomson scattering generates linear 
polarization with a polarization tensor j.'ill I32| 

"Pab(e) = -^drV <0 V 6> e 2 , (17) 

where dr is the differential optical depth, and the covariant derivatives are 
in the surface e 2 = 1. Polarization is only generated from the quaudrupole 
component of the incident radiation, and can be seen to be purely electric 
quadrupole in character. If the incident radiation is polarized, there is an addi- 
tional contribution to the polarization produced due to the electric quadrupole 
of the incident polarization (0 2 should be replaced by 2 — 12Pe j2 [2JC21)- 
For the CMB, the linear polarization we observe (at x = 0) along a line 
of sight ft is that produced at last scattering along the radiation direction 
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e = — n, but at a direction-dependent comoving position x = Arjfi. (We are 
assuming the Universe is flat for simplicity, and Ar) is the conformal time 
between last scattering and the present.) The character of the observed po- 
larization thus depends on the spatial distribution of the polarization at last 
scattering which, in turn, depends on the type of perturbations (e.g. density 
or gravitational waves) we consider. 

3.1 Density Perturbations 

If we consider a single plane- wave density perturbation of comoving wavenum- 
ber k, the peculiar velocity of the electron-baryon plasma around recom- 
bination is in the direction of the wavevector k. Over the mean free time since 
the previous scattering, a temperature quadrupolc 

2 (e) ~ ^pTcMBV^t^eV cx T CMB £ p kP 2 (k • e) (18) 

is produced to leading order in £ p k where £ p is the mean free path. The 
quadrupole is azimuthally symmetric about the wavevector, and Thomson 
scatters to generate linear polarization parallel to the projection of k on the 
sky (see Fig. . The polarization tensor at last scattering has a spatial de- 
pendence exp(ifc • x) appropriate to the plane wave, and this leads to a further 
modulation of the polarization observed today by exp(iAryfc • fi). This mod- 
ulation transfers polarization from the / = 2 to higher multipoles (with most 
appreciable power appearing at I = fcAry), but preserves the electric charac- 
ter of the polarization [111 ITS] . To see that this is reasonable, note that the 
modulated polarization field has its polarization direction either parallel or 
perpendicular to the direction in which the polarization amplitude is chang- 
ing PSj- Comparison with Fig. shows that such a pattern is pure electric 
polarization. This important observation, that linear density perturbations do 
not produce magnetic polarization, can also be understood from the fact that 
the spatial distribution of V^w bj ) is curl-free [3*2| . 

An example of the polarization power spectra produced by adiabatic den- 
sity perturbations is plotted in Fig. They were computed with the Boltz- 
mann code CAMB 34]. The -E-mode power peaks around I ~ 1000, corre- 
sponding to the angle subtended by the width of the visibility function at 
recombination. On larger scales the polarization probes the electron-baryon 
velocity at last scattering, as described above. Acoustic oscillations of the 
plasma prior to last scattering imprint an oscillatory structure on the an- 
gular power spectra of the CMB observables. Fourier modes of the photon 
density that are caught at the extrema of their oscillation at last scatter- 
ing produce the peaks in the temperature-anisotropy power spectrum, but 
troughs in the polarization since then vanishes for such modes. Modes 
caught at cither the extrema or the midpoint of their oscillation give zeroes in 
the temperature-polarization cross-correlation. Large-angle polarization from 
the last scattering surface is very small: the optimal wavelength for generating 
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Fig. 4. Generation of linear polarization by Thomson scattering for a single Fourier 
mode of the density perturbations, following the treatment of The temperature 
anisotropics are azimuthally symmetric about the wavevector; the quadrupole part 
is shown on the left with cold lobes along the indirections. The linear polarization 
that is produced is then a pure I — 2, m = electric polarization pattern locally 
(right). To get the observed polarization after free streaming from last scattering, 
the local polarization is modulated by the plane wave exp(ifc • nAij). 



power at scale I is l/Arj, but this is too small to allow significant generation of 
quadrupole temperature anisotropies over a mean free path before last scat- 
tering. (Large-angle polarization is further suppressed since V\, is small on 
super-Hubble scales.) The increase in polarization on large scales in Fig. [2] 
is due to reionization [TJ3 . Re-scattering of the temperature quadrupole gen- 
erates linear polarization peaking at multipolcs roughly twice the ratio of 
the conformal look-back time to reionization to the conformal time elapsed 
between recombination and reionization. The amplitude of the polarization 
produced at reionization is proportional to the number of photons that scat- 
ter (i.e. the optical depth to reionization, taken to be r = 0.15 in Fig.^J. 

Although linear density perturbations do not produce £?-mode polariza- 
tion, this is produced by a number of second-order effects. The most notable 
is weak gravitational lensing by large-scale structure jHS], which effectively 
re-maps the polarization on the sky according to the lensing deflection field. 
This can distort a curl-free pattern into one with curl and so generate mag- 
netic polarization. The power spectrum of this effect is also included in Fig. [21 
On large scales the spectrum is almost white with Cf ~ (1.3 nK) 2 . The lens- 
induced B-mode spectrum peaks at I ~ 1000 corresponding to the peak in 



3.2 Gravitational Waves 

A stochastic background of gravitational waves, as predicted in all infla- 
tionary models, can also influence the CMB. Locally, a gravitational wave 
causes a quadrupole anisotropy in the expansion of space, so photons arriv- 
ing at a scatterer from different lines of sight will suffer different redshifts 
since their previous scattering event. The longest wavelength gravitational 
waves are geometrically most efficient at generating a quadrupole moment. 
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For wavelengths long compared to the mean free path, the quadrupole mo- 
ment 62(e) oc —IpTcMBKj^^/^, where h^ is the transverse, trace-free per- 
turbation to the Robertson-Walker metric, 

ds 2 = a 2 (r])[dr] 2 - (<% + h , , )<\.r' a.r ! , (19) 

and an overdot denotes differentiation with respect to proper time. However, 
the shear of the gravitational wave hij is zero for wavelengths well outside the 
horizon, and undergoes damped oscillation inside the horizon, ft follows that 
the temperature quadrupole at last scattering is dominated by modes with 
wavelength of the order of the horizon size, and this sets the characteristic scale 
(I ~ 100) on which the polarization power spectra peak (see Fig.EJ. Thomson 
scattering of the temperature quadrupole from gravitational waves generates 
a polarization tensor at last scattering V a b oc [/i b] TT , i-e. proportional to 
the trace- free projection of the shear onto the sphere e 2 = 1. For a single 
plane wave with wavevector along the z-axis, the quadrupole moment and 
local polarization from the plane wave have only m = ±2 modes (see Fig. \5\ 
and |31| for further details). The observed polarization from a plane wave is 
further modulated by exp(iA?]fc • n) in free-streaming from last scattering. 
Modulating the polarization in Fig. [S] gives a polarization pattern that locally 
looks electric along the lobes of the temperature quadrupole, and magnetic 
in between jSSJ- It follows that both electric and magnetic polarization are 
generated by gravitational waves |lllll2j . and with roughly equal powers |31| . 
The generation of B-mode polarization from gravitational waves can also be 
understood as a consequence of the non-vanishing (spatial) curl of the shear 
hij at last scattering |3~2*] . 

The power spectra of the CMB observables from a scale-invariant back- 
ground of gravitational waves is shown in Fig. [2 As for density perturbations, 
large-angle polarization is generated by reionization. The amplitude of any 
gravitational wave background is currently only weakly constrained by CMB 
observations. A recent analysis, combining CMB observations and data from 
the Sloan Digital Sky Survey (SDSS), constrained the tensor-to-scalar ratio 
r — a measure of the ratio of the primordial power in gravitational waves to 
curvature (density) perturbations - to r < 0.36 at 95% confidence in a seven- 
parameter flat model. This is similar to the value adopted in Fig. [3 and shows 
that the gravitational wave contribution to the temperature anisotropies and 
electric polarization is sub-dominant. 

4 Current Status of CMB Polarization Measurements 

The first detection of polarization of the CMB was announced in September 
2002 |SJ. The measurements were made with DASI, a compact interferomet- 
ric array operating at 30 GHz, deployed at the South Pole. The DASI team 
have now analysed three years of data; they report a detection of -E-mode po- 
larization at high significance (6.3er) and at a level perfectly consistent with 
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Fig. 5. As in Fig.^Jbut for a Fourier mode of the gravitational wave perturbations. 
The quadrupole part of the temperature anisotropy now has m = ±2 modes (left); 
the hot lobes correspond to the directions in which space is locally being contracted 
by the wave. The linear polarization produced by Thomson scattering (right) is again 
electric and quadrupole in character at the last scattering surface. Unlike density 
perturbations, magnetic polarization is now generated by free streaming since the 
polarization field modulated by the plane wave has a non-zero curl. 

predictions based on temperature-anisotropy data |36| . The only other re- 
ported direct detections to date are from the Cosmic Anisotropy Polarization 
MAPper (CAPMAP; 0), a heterodyne correlation polarimeter operating at 
90 GHz from New Jersey, and from the Cosmic Background Imager (CBI; [§]), 
an interferometer operating in the range 26-36 GHz from Chile. The results 
on Cf from maximum-likelihood power spectrum analyses by the DASI, CBI 
and CAPMAP teams are shown in Fig. along with the prediction based 
on the best-fit to all the first-year WMAP data. Although the data are very 
noisy and broad-band (due to the limited sky coverage), the agreement with 
the prediction is striking. In addition, the TE cross-correlation has been mea- 
sured with good precision by WMAP 0, and has also been detected by DASI 
(see Fig. 01 . 

The -B-mode polarization of the CMB is currently subject only to weak 
upper limits. The best limit on large scales, I ~ 10, is from the POLAR exper- 
iment , while the best limits on smaller scales are from DASI [35]. These 
(95%) upper limits are shown in Fig. along with the predicted .B-mode 
spectrum in a model with tensor-to-scalar ratio r = 0.38, close to the current 
95% limit from temperature anisotropies and SDSS |38| . Current direct con- 
straints on the amplitude of gravitational waves from i?-mode data are clearly 
much weaker than this. Furthermore, upper limits on the B-mode power on 
small scales are at least two orders of magnitude above the level expected 
from gravitational lcnsing. 
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Fig. 6. Current polarization detections and best upper limits (as of January 2005). 
The top plot shows Cj B measurements (top) and Cf (bottom). The points with 
errors are from the first one-year data release from WMAP. The error boxes 
are the flat band-power results from DASI, CBI and CAPMAP data, centred on 
the maximum-likelihood band power and spanning the 68-per cent intervals. The 
solid lines are the predicted power from the best-fit model to all the WMAP data. 
The bottom plot shows the upper limits on Cf from DASI, CBI, CAPMAP and 
POLAR. The solid line assumes a tensor-to-scalar ratio r = 0.38. 

5 What do we learn from CMB polarization? 

Electric polarization primarily traces the velocity of the electron-baryon 
plasma at recombination due to density perturbations. It thus provides com- 
plementary information to the temperature anisotropies which are sourced 
mainly by the photon density variations at last scattering on sub-degree scales. 
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Polarization observations can thus be expected to bring significant improve- 
ments on some of the 'acoustic' cosmological parameters [5], such as the phys- 
ical density of baryons and cold dark matter, although current polarization 
data is not yet of sufficient quality relative to the temperature data to see this 
effect 6 . Furthermore, the relative phase of the acoustic peaks in the temper- 
ature and electric-polarization power spectra (and the cross-correlation) can 
be used for a largely model-independent test of the physics of acoustic oscil- 
lations. For both adiabatic and isocurvature perturbations the acoustic peaks 
in Cj and Cf should be in antiphase. (Admixtures of these perturbations 
can, however, partially destroy this coherence.) A recent analysis by the CBI 
team found the phase of the oscillations in the E'-mode spectrum to be 
fully consistent with expectations for adiabatic models. The spectrum Cf E 
from the one-year WMAP data [H] has also been used to improve constraints 
on the amplitude of any isocurvature fluctuations, e.g. HU1 HT] . 

Polarization on scales of tens of degrees (I < 10) can be used to probe 
reionization. Significantly, it breaks the degeneracy between the amplitude of 
density perturbations, A,, and the optical depth to reionization, r, that leaves 
only A s e~ 2r well-determined from the temperature anisotropics. (The power 
spectra of the CMB observables are reduced by e _2r on scales corresponding 
to Fourier modes that are sub-horizon at reionization.) The one-year WMAP 
measurement of Cf E shows a significant excess power on large scales that is 
attributed to early reionization (11 < z rc < 30) giving a significant optical 
depth r = 0.17±0.04 (68% interval) [Hj- The amplitude of curvature (density) 
perturbations (at 0.05 Mpc) is then constrained to A s = (2.7±0.3) x 10 -9 from 
the CMB alone [3U] ■ Looking further ahead, the fine details of the large-angle 
polarization power can, in principle, distinguish different ionization histories 
with the same optical depth, although this is hampered by the large cosmic 
variance at low I 021 E3| • 

Perhaps the greatest ambition for CMB polarization observations is the 
detection of the -B-mode signal from gravitational waves. Although gravita- 
tional waves do generate temperature anisotropics and electric polarization, 
they do so only on large scales where the cosmic variance of the dominant 
density perturbations is large. Indeed, a perfect temperature-anisotropy ex- 
periment can only detect r > 0.07, even assuming all other parameters are 
known. Supplementing the temperature anisotropics with electric polarization 
makes only a modest further improvement: r > 0.02. Although the current 
upper limit on r from B-mode polarization is much weaker than that inferred 
from the CMB temperature anisotropies (see Fig. the former route ulti- 
mately has the potential to probe much lower values as thermal noise levels 
improve. Assuming that astrophysical foregrounds can be removed adequately, 
the ultimate limit to a detection of the gravitational wave signal is set by the 
B- modes generated by gravitational lensing ^] • 

A detection of a stochastic background of gravitational waves would have 
significant implications for fundamental physics. In inflation models, the power 
produced in gravitational waves as a function of wavelength is a direct mea- 
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sure of the expansion rate of the Universe at the time that wavelength was 
stretched beyond the horizon. In slow-roll inflation in a potential V, the power 
in gravitational waves thus constrains the 'energy scale of inflation', V 1 ^ 4 , via 

V 1/A = 3.33 x 10 16 r 1/4 (yl s /2.3 x 1(T 9 ) 1/4 GeV . (20) 

The current 95% limit r < 0.36 |2Hj constrains V 1 / 4 < 2.3 x 10 16 GeV. By way 
of contrast, the recently-introduced cyclic model |47j would produce negligible 
gravitational waves as the fluctuations are produced during a phase of very 
slow contraction 0H| • A detection of cosmological gravitational waves would 
effectively rule out such models. In addition to constraining the energy scale of 
inflation, measuring r is important for constraining the shape of the inflaton 
potential. For the simplest models of slow-roll inflation, the observables r and 
n s (the spectral index of curvature perturbations) are both required to specify 
the first two derivatives of the inflaton potential. Without the tight limits on 
r from future -B-mode searches, constraints on inflation models will be largely 
degenerate along the direction of constant n s . The most recent constraints on 
r and n s from the CMB and galaxy surveys already rule out some popular 
models such as V oc <fi A |38| . 

There are several other potential sources of .B-mode polarization that could 
dominate over any signal of primordial gravitational waves from inflation. A 
sub-dominant background of cosmic strings - currently enjoying a resurgence 
of interest following the realisation that they are naturally produced in some 
models of brane inflation |49j - would provide a seed for gravitational waves 
(and vortical, or vector, modes) on the last scattering surface. Both give rise to 
i3-mode polarization, and the resulting power spectrum is expected to peak 
on sub-degree scales [SOI. It has been argued that a statistically-isotropic, 
but highly non-Gaussian, pattern of .B-mode polarization with power peak- 
ing on small scales would be a smoking-gun signature for a background of 
cosmic strings |51) . In addition, primordial magnetic fields have several in- 
teresting effects on the CMB and its polarization, although the theory is still 
rather under-developed. The anisotropic stress of a stochastic field also excites 
gravitational waves and vector modes (e.g. [52]) leading to (non-Gaussian) B- 
mode polarization. Faraday rotation by either large-scale coherent j53J or 
stochastic fields [HE] can also rotate electric polarization into magnetic. The 
hope is that these effects could be separated on the basis of their non-Gaussian 
statistical properties, and, in the case of Faraday rotation, by its frequency 
dependence. 

Finally, we note that weak gravitational lensing effects on the CMB are a 
valuable source of additional cosmological information, most notably on the 
properties of dark energy and neutrino masses [HE]- As we describe further 
in Sect. 16.61 CMB polarization is particularly useful here since intrinsically it 
has more power on small scales than the temperature and so is more sensitive 
to the small-scale power in the lensing deflection. 
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6 Polarization data analysis 

The CMB polarization signal is small. Current measurements are consistent 
with the expected r.m.s. signal ~ 6 fiK, which should be compared with the 
~ 120 /iK r.m.s. of the temperature anisotropies. The low signal level requires 
high thermal sensitivity and demands exquisite control of systematic effects. 
The impact on the final science products of any systematics that cannot be 
designed out of the instrument must be carefully assessed (usually with a 
large number of simulations), and any that are found to be significant must 
be accounted for in the data analysis. As noted in the Introduction, this in- 
evitably makes the analysis of polarization data more instrument-specific than 
for temperature- anisotropy experiments. Despite this, there are a number of 
common steps in the analysis pipeline of any polarization experiment: map- 
making from time-series data; astrophysical foreground removal; E-B mode 
separation; power spectrum estimation; and parameter estimation. (Separa- 
tion of E and B modes can also be performed statistically during the esti- 
mation of power spectra, and this has been the case for all analyses of real 
data to date.) In the following subsections we briefly summarise most of these 
key steps, but focus on E-B separation and power spectrum estimation in 
particular. While the other steps are, of course, crucially important, the anal- 
ysis methods employed there are relatively straightforward generalisations of 
those developed for the temperature anisotropies, for which excellent reviews 
already exist (see e.g. |S7])- 

6.1 Modelling the polarimeter response 

We begin by looking at how to model the response of a polarimeter, i.e. the 
relation between the observed data in the time domain to the fields on the 
sky. The most general (instantaneous) linear response of a polarimeter to the 
polarized sky brightness at each frequency is of the form 

s oc / (77 + QQ + UU- VV) dn (21) 

Here, I(ii 1 v), Q(n,i/), U{n,v) and V{n,v) are the Stokes brightness pa- 
rameters at frequency v, and the equivalent quantities with tildes refer to the 
instrument response. The latter parameters are determined by the instrument 
optics and receiver. They inherit a time-dependence from the scanning of the 
instrument on the sky - 7(n, v) and V(n, v) rotate with the instrument as 
scalar functions and Q(n, v) and U{h, v) as components of a spin-2 function 
- but also from any modulation scheme employed to ameliorate the effects of 
low- frequency noise in the instrument. The signal contribution to the data is 
given by further integrating i|21|) over the spectral response of the instrument. 

A useful idealisation is an instrument with ideal optics and receiver, and 
a narrow spectral bandpass. For a dual-polarization system, such as a pair 



18 Anthony Challinor 



of polarization-sensitive bolometers on the high-frequency instrument of the 
Planck satellite, the beam patterns on the sky for the two polarization modes 
should have field directions that are orthogonal and constant, and the am- 
plitude of the response should be azimuthally symmetric. In the case of the 
bolometric system, the ideal instrument response is then fully determined by 
I(n), which is a function of 9 only. The response to circular polarization van- 
ishes and Q ± iU = —I(9)e ±21 ^ for the y-polarization, and minus this for 
the ^-polarization. If we now rotate the instrument so the optic axis is along 
the direction n, and the local y-polarization direction is at an angle ip to the 
0-dircction there, the signal received in the y-polarization simplifies to |58| 

J (II + QQ + UU- VV) drc = J 8 (n) - Q s (n) cos 2^ + U s (n) sin 2ip , (22) 

where the beam-smoothed fields are 

I s (n) =J2WiT lm Y lm (n) , (23) 
im 

(Q s ± iU s )(h) = ME T iB) lmT2 Y lrn (n) . (24) 

l in 

Here, the Wi is proportional to the m = multipoles of 1(9) and 2W1 to the 
m = 2 electric multipoles of Q±iU. For a Gaussian beam of dispersion <j <C 1, 
we have beam functions 

WlKe -l{l+X)^/2 and ^ M e _ [Km) -4] CT =/2 . (25 ) 

For I 3> 1 the polarized and total intensity beam functions are very nearly 
equal. The form of (|22|l is equivalent to a point sampling of the smoothed 
fields along the optic axis with a polarizer in orientation ip. The signal in the 
a;-polarization is got by replacing ip by ip — tt/2. The generalisation of 11' L'l) to 
other ideal polarimeters is straightforward: in each case the signal remains a 
linear combination of the beam-smoothed fields. 

6.2 Polarized map making 

The aim of map-making is two-fold: (i) to produce pixelised images of the 
temperature and polarization fields on the sky from the observed data; and (ii) 
to compress the time-stream data to a more manageable size, while preserving 
as much of the cosmological information present as possible. We shall only 
consider the ideal instrument model described above; for a case study of some 
of the real-world problems that must be solved in the map-making process for 
total-intensity data see jSH]. For the ideal model, the signal contribution to 
the data is linear in the beam-smoothed fields. Instrument noise is assumed 
additive, in which case our simple model for the data from the ith detector 
at time t becomes 
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4 = 4 p (I SjP - Q s , p cos 2^t + U s , p sin 2$) + n\ . (26) 

Here, A\ p is the detector pointing matrix and is unity if the beam centre of 
detector i points towards pixel p at time t, and is zero otherwise. The ith 
detector at time t is sensitive to polarization along a direction at angle ipl 
to the ^-direction on the sky, and has instrument noise n\. We assume that 
all detectors have the same narrow spectral bandpass, and the same angular 
resolution, so that we can combine the signals of all detectors into a single 
vector d and the instrument noise into n. We attempt to reconstruct the 
pixelised beam-smoothed fields, which we represent by the vector s, by a 
regularised inversion of d = As + n, where the (it)p element of A is A l tp . 

If the noise is Gaussian, with known covariance (nn T ) = N, the maximum- 
likelihood solution, 

s = (A T N _1 A) _1 A T N _1 d , (27) 

is optimal in the sense of being the minimum-variance unbiased estimate. 
The covariance of the errors is ((s — s)(s — s) T ) = (A T N _1 A) _1 . If the 
CMB fields are Gaussian, the maximum-likelihood map contains all of the 
cosmological information present in the time-stream data. In the language of 
statistics, the maximum-likelihood map is a sufficient statistic for parameter 
estimation, i.e. the data d only enters the posterior probability of the cho- 
sen (cosmological) parameters a. given the data, Pr(a:|d), through s. Forming 
the maximum-likelihood map can thus also be viewed as lossless compres- 
sion of the time-stream data. The question of the feasibility of computing the 
maximum-likelihood map is discussed in discussed in detail elsewhere in this 
volume by Borrill. Here we make just a few general remarks. 

• Any matrix can be used in place of the noise covariance matrix N in 1271) 
and the solution is still unbiased. However, if the matrix used does not 
properly capture the noise correlations described by off-diagonal elements 
of N, the resulting maps will have correlations along the scan directions 
('stripes'). 

• Brute-force evaluation of s is prohibitive for large datasets so approximate 
methods must be used. For stationary noise, N is circulant and can be 
applied efficiently in the Fourier domain. Furthermore, direct inversion of 
A T N _1 A can be avoided by solving A T N _1 As = A T N _1 d with iterative 
techniques. 

• For stationary noise, the power spectrum (the diagonal elements of N in 
Fourier space) can be estimated from the data if it is not known a priori. 
Errors in the noise estimation increase the errors in the map. 

• Maximum-likelihood map-making is easily generalised to include template 
fitting and marginalisation. For example, the unknown amplitude of some 
additive systematic effect whose 'shape' is known can be marginalised over 
during the map-making stage. 

• Fast, instrument-specific alternatives to maximum-likelihood map-making 
have been proposed, e.g. the 'destriping' technique for Planck |6U| . 



20 Anthony Challinor 



6.3 Astrophysical foreground removal 

The dominant diffuse polarized foregrounds in the frequency range 20- 
400 GHz are expected to be Galactic synchrotron and vibrational (thermal) 
dust emission. Synchrotron radiation in a uniform magnetic field is linearly 
polarized. The polarization fraction depends on the energy spectrum of the 
electrons, but can plausibly be as high as 70%. Thermal dust emission will 
also be linearly polarized if non-spherical grains are aligned in a magnetic 
field. Any variation of the Galactic magnetic field along the line of sight will 
tend to depolarize the radiation reducing the polarization fraction. This depo- 
larization tends to cancel the increased emission along the plane making the 
polarization dependence on Galactic latitude weaker than in total intensity. 
Synchrotron and dust emission can be expected to produce both electric and 
magnetic polarization in similar amounts |61j . 

Synchrotron radiation is expected to be the dominant polarized foreground 
at frequencies below ~ 100 GHz. At the time of writing, the only well-surveyed 
part of the sky for which polarized data is available is within the Galactic 
plane, and then only at frequencies below 2.7 GHz (e.g. The only high- 

latitude observations available are very patchy and most are not useful 
for assessing foreground contamination in clean regions of the sky that will 
likely be the target of future deep polarization observations. (But see [HI] for 
a useful exception.) Analysis of the 2.7 GHz data of the Southern Galactic 
plane from gives E- and i?-mode power spectra going as I -1 - 5 [55], but 
care should be taken in extrapolating this to higher latitudes. A power law 
index of —1.5 is shallower than the large-angle polarization for both density 
perturbations and gravitational waves. Modelling at higher latitudes, based 
on total intensity observations (e.g. [US]), suggests that at frequencies around 
100 GHz, synchrotron emission will not be a major foreground for £ , -modc 
polarization, but that its removal will be essential for large-angle S-modc 
searches. We eagerly await the release of polarization maps from the two-year 
WMAP data; these should greatly improve our knowledge of high-latitude 
synchrotron polarization. 

The power spectra from diffuse Galactic polarized dust emission have been 
measured recently at 353 GHz by the Archeops team [HH] over 20% of the sky. 
Extrapolating their measurements to 100 GHz, they find for 3 < I < 70 that 
(I + l)Cf' B /2ir < 0.2 /iK 2 , excluding data with 5° of the Galactic plane. 
Removal of the dust contamination is thus also likely to be necessary for 
future large-angle B-mode searches. 

A variety of methods have been developed to remove diffuse foreground 
emission from CMB maps using multi-frequency data; see the contribution by 
Delabrouille in this volume for a review. Many of these were developed for total 
intensity data (temperature), but they generalise to polarization straightfor- 
wardly when applied directly to Stokes maps. Of these, some, such as Weiner 
filtering |67) . require knowledge of the frequency spectra of all the compo- 
nents present, and their angular power spectra. Others attempt to derive this 
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information from the data by assuming statistical independence of the com- 
ponents Finally, linear techniques, based on template fitting P, or 
internal spectral combinations that preserve the CMB signal while minimising 
the r.m.s. of the combined field, can also be applied very effectively to polar- 
ization data provided that the spectral properties of the contaminants do not 
vary too greatly over the survey area. The results of applying these techniques 
to simulated polarization data are very encouraging (e.g. |7()j). but we must 
await further observations to assess the true extent to which foregrounds will 
compromise the scientific returns of CMB polarization surveys. 

Extra-galactic radio sources are a significant polarized contaminant at 
small angular scales. They contribute equally to E- and B-mode power and, 
for Poisson-distributed sources, have white angular power spectra. An effec- 
tive way to deal with these is simply to exclude pixels within an observing 
beam of known sources above some flux cut. Since the degree of polarization is 
typically rather low (< 10%), source contamination for iJ-mode polarization is 
relatively less troublesome than for temperature anisotropies. Having Stokes 
maps with a large number of excluded pixels does complicate the problem of 
separating E and B modes, and this is worse for low-resolution observations 
since the fraction of the area that is removed is then larger. 

6.4 Map-level E and JB-mode separation 

Separating the recovered polarization field into its electric and magnetic parts 
amounts to solving for the potentials Pe and Pb in ©. This can be done 
trivially if the observation covers the full sky by using the orthogonality of 
the electric and magnetic tensor harmonics to recover the multipoles Ei m 
and E>i m as in Q. However, problems arise for observations over only part 
of the sky since the decomposition V a b — V^V^Pe + e c ^ a V & )V c Pe is then 
not unique [721 171). Even if fully-sky observations were available, we would 
want to remove certain regions (e.g. the Galactic plane) from the subsequent 
analysis, and this is most naturally done before the non-local decomposition 
into E and B modes. 

To see how the ambiguity comes about, consider attempting to solve for 
Pe and Pb by recalling that the electric contribution to the divergence V b "P a fc 
is curl-free and the magnetic part is divergence free, i.e. 

2V a W b V ab = V 2 (V 2 + 2)P E , (28) 
-2e ac V c S7 b V ab = V 2 (V 2 + 2)P B . (29) 

Solving these equations over part of the sphere requires one to specify bound- 
ary data, namely the value of the potentials and their normal derivatives on 
the boundary [JT|. But this information is not known since it is the poten- 
tials that we are attempting to solve for. We can, of course, impose boundary 
conditions by force, in which case we should ensure that we do not produce 
spurious E or B modes, i.e. we should recover Pb = if there are no B- 
modes present. This will be true if we solve (|28f) and (|29|) with the boundary 
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conditions that both the potentials and their gradients vanish on the bound- 
ary |72M71| . In this manner, we obtain a unique solution for P E and P E , but 
they are systematically wrong by AP^ and APg which are the unique solu- 
tions of V 2 (V 2 + 2) APe = with AP E = P E and V a AP E = V a P E on the 
boundary (similar equations hold for APe). Reconstructing the polarization 
field from the recovered potentials, we see that this differs from the original 
field by AV a b = V( a Vb)AP_E + e c ( a Vf,)V c APB. This part of the polarization 
field, which satisfies V V b P a & = and e ac V c V b AP a (, = 0, constitutes the 
ambiguous part since it cannot be classified in terms of electric and magnetic 
modes. Some examples of these ambiguous modes are illustrated in |71l I73| . 
A useful vector analogy to keep in mind is that the ambiguous modes are like 
the electromagnetic fields of (static) charges and currents that are localised 
outside the observed region. 

By excluding the ambiguous modes from subsequent analysis, we are 
throwing away some information. As we might expect, the ambiguous modes 
tend to be localised near the boundary, and a rough accounting of the frac- 
tional loss of modes as a function of scale is given by the ratio of the product 
of the boundary length and the coherence length to the observed area. 

Practical methods for performing the separation described above are given 
in [121 El EH- The basic idea is to avoid differentiating the data by instead 
constructing tensor bases that can be used to project out the unambiguous 
E and B modes. Consider, for example, isolating magnetic polarization. If 
we construct a symmetric, trace-free tensor W Ea b = V / 2e c ( a ^h) V C M^ from a 
scalar field W that vanishes, along with its gradient, on the boundary but is 
otherwise arbitrary. Contracting with the polarization tensor and integrating 
by parts gives [72] 

/ W^Vab dn = -2 y W*e b c V c V a V ab dn , (30) 

where there are no boundary terms on account of the conditions imposed on 
W. It follows that the integral over the observed area is a non-local measure 
of the P-mode component. Similarly, the tensor W Ea b = v / 2V( a V b )W, with 
the same scalar field W, extracts a measure of the P-mode component. The 
construction of a complete set of such tensors can be performed in spherical- 
harmonic space [731731 (see also Sect. 16 .5|) . or directly in pixel space [7T]. An 
orthonormal tensor basis can be constructed by taking the W to be eigenfunc- 
tions of V 2 (V 2 + 2) (with the boundary conditions given above). Eigenfunc- 
tions with different eigenvalues are then orthogonal over the observed area, 
and the same is true of the tensors derived from them. Examples of W E ab 
basis functions for a circular patch of the sky are shown in Fig.[7J In this case 
the tensor bases can be constructed with W an eigenfunction of d/d^. For 
m > 2 exactly two P-modes are lost per to; for m = only one is lost; and 
for rn = none are lost. The basis functions in the figure are constructed to 
maximise the expected power in that mode. 
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Fig. 7. Magnetic tensor basis functions WBab for projecting onto pure magnetic- 
polarization modes for an azimuthally-symmetric sky patch of radius 10° 1721 . Only 
m — 0, 1 and 2 modes are shown. In each case only the three modes with the largest 
expected power are shown (with the left-most having greatest power). 

The loss of information due to the ambiguous modes has interesting im- 
plications for survey design. For a fixed integration time, the 'best' strategy 
for measuring the variance of a Gaussian, statistically-isotropic field, in the 
limit of low signal-to-noise, is to survey only a small field since the error on 
the variance scales as the square-root of the survey area. (Of course, this is 
not a wise strategy as you may be unlucky and choose a field where the sig- 
nal happens to be very low.) Applying this to searches for gravitational waves 
with B-mode polarization, the implication is that a small survey is best. How- 
ever, the coherence length of the desired signal is around a degree so the loss 
of information to ambiguous modes becomes significant for small fields. The 
optimum field size on this basis turns out to be a circular patch with radius 
around 10° [72] • This argument ignores the effect of weak gravitational lensing. 
If we include the lens-induced B modes as an additional Gaussian noise term, 2 



2 The B modes produced by lensing are not Gaussian. Their (connected) four- 
point function does have a significant effect on the variance and covariance of the 
estimated B-mode power spectrum |74j . 
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i.e. we make no attempt to clean out the lensing contribution as described in 
Sect. 16.61 the optimum size is increased further such that the contributions to 
the error on the variance from thermal noise and lensing sample variance are 
roughly equal. 

We end this subsection by noting that coherent E-B separation is not 
strictly necessary for constraining cosmology. A statistical separation can be 
attempted during power spectrum estimation (see Sect. l6.f)|) . and for the opti- 
mal maximum-likelihood method this should be very effective. However, there 
are good reasons to perform a coherent, map-level separation as well. For ex- 
ample, visualising the B modes will likely be a useful diagnostic of unknown 
systematic effects and foreground residuals. 

6.5 Power spectrum estimation 

The statistical properties of Gaussian random fields are fully specified by their 
power spectra. For this reason, power spectrum estimation is an important 
step in the analysis of CMB data. Ideally, this involves mapping out the like- 
lihood Pi(s\Cp-), where s are the observed CMB fields, as a function of the 
power spectra C/* (R = T, E, B or TE). In practice, what is normally done 
is to approximate the likelihood as a Gaussian function of the power spec- 
tra, in which case all of the information contained in the original CMB maps 
is distilled into the position of maximum likelihood, and the curvature of the 
likelihood there 75 . In this process, we achieve a very significant compression 
of the CMB maps into, at most, a few thousand maximum-likelihood C;S, and 
their associated covariances. Even this programme is not achievable exactly for 
large data sets, so a number of faster techniques based on estimators have been 
suggested. Here, the idea is to construct (unbiased) estimators of the power 
spectra, and to use their covariances - or, better, their full sampling distribu- 
tions _ to constrain cosmology. A particularly fast class of estimators are 
the heuristically- weighted, quadratic, or 'pseudo-CV estimators (e.g. [771 178] ). 
For an interesting critique of pseudo-C; and maximum-likelihood methods for 
temperature anisotropies, see [7JJ]. For polarization, both methods can be ap- 
plied directly to maps of the Stokes parameters, in which case one must rely 
on a statistical separation of the E and B modes, or after map-level E-B 
separation. 

An important exception that demands a treatment beyond the power spec- 
trum is the effect of weak gravitational lensing. Special techniques have been 
developed for this problem which allow reconstruction of the underlying Gaus- 
sian CMB and deflection fields (see Sect. I6.6fl . The power spectra of these 
reconstructed fields then contain all of the statistical (hence cosmological) 
information present in the non-Gaussian, lensed fields. 

Maximum-likelihood methods 

Maximum-likelihood methods, and the problems one faces in their practical 
implementation, are reviewed extensively elsewhere in this volume by Bor- 
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rill. Here, we briefly summarise the technique. Maximum-likelihood methods 
have been used extensively in the analysis of real polarization data, e.g. from 
DASI and CBI 0. 

For the direct analysis of Stokes maps (i.e. assuming no E-B separation 
has been performed) we combine the observed Stokes parameters in a vector 
s = (Ip,Q p ,Up) T . The covariance matrix of this vector, C = S + N, has 
contributions S from the signal and N from the noise. The former is given 
essentially by the temperature and polarization correlation functions <|12|1 . 
but with the following modifications: (i) the power spectra are multiplied by 
the appropriate product of the beam functions Wi and 2W1 to account for 
convolution with the instrument beam; and (ii) there are additional factors to 
account for the rotation of the Stokes parameters from the polar basis to the 
geodesic basis defined by the points p and p 1 . The noise contribution N includes 
the projection of the instrument noise onto the map, which will generally be 
modified by the foreground removal process, and an estimate of the covariance 
of any residual foreground errors. Both S and N have (3iVpi x ) 2 elements; S is 
dense in the pixel domain but is sparse in spherical-harmonic space, while N is 
typically sparse in the pixel domain. If map-level E-B separation has already 
been performed, the data s can be taken as the unambiguous E and B modes, 
and these can be analysed independently (although the noise may correlate the 
errors in the estimated E- and -B-mode power spectra). Assuming Gaussian 
signals and errors, the likelihood function is 



Pr(s|C/ 



1 



v /det(2 7 rC) 



exp(-s T C" 1 s/2) . 



(31) 



Maximising the likelihood is equivalent to maximising the posterior probabil- 
ity Pr(C i R |s) if a uniform prior on the power spectra, Cf, is adopted. 
We now make some general remarks on this procedure. 

• For full-sky coverage, and no noise, the log-likelihood is (up to an irrelevant 
constant) 



21nPr(s|Cf) = £>i + 1) 



E 



Er<T 



c?c, 



2Cf E Cf E ) 



■ In A ; + 



Cf 

cf 



\^Cf 



(32) 



where A; = CfCf — (Cf E ) 2 . The likelihood is approximately a Gaussian 
function of the Cf at high I. The maximum-likelihood point coincides with 
the simple estimators e.g. Cf = J^m \Eim\ 2 /(2l + 1) at low I. 
Generally, the maximum of the likelihood can be located by solving 
d\nPr(s\Cf)/dCf = 0, with e.g. the Newton-Raphson method. The data 
enters the gradient of the log-likelihood, and hence the maximum likeli- 
hood point, in the form s T C^ 1 (dS / dCf jC^ 1 s . This is proportional to the 
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simple estimate Cf 1 , but with the spherical harmonic transforms taken 
only over the observed region, and with the data first weighted non-locally 
by C '. 

• Approximating the likelihood as Gaussian requires the computation of the 
curvature of the likelihood. At any Cf 1 , the expectation value of the second 
derivatives of (minus) the log-likelihood there, under the assumption that 
these are the true power spectra, defines the Fisher matrix. The inverse 
of the Fisher matrix, evaluated at a smoothed version of the maximum 
likelihood power spectra, is often quoted as the error on the power spectra. 
For the ideal case, where the log-likelihood is given by l|32[) . the inverse of 
the Fisher matrix has elements given by the cosmic variance formulae of 
Sect.O 

• A brute-force implementation to locate the maximum likelihood, and eval- 
uate the Fisher matrix, is not feasible for mega-pixel data sets. Both the 
operation count, 0(iVp ix ), and the storage requirements are prohibitive. 

• Additional regularisation is required for small surveys since the problem is 
then degenerate if we attempt to solve for every multipole A common fix 
is to solve for bandpower amplitudes |75| . For narrow bands, approximat- 
ing the 1(1 + l)C;s as piecewise constant is common. As the polarization 
power spectra are steeper functions of I than C ; T , particularly on large 
scales, shaped bandpowers are often employed in polarization analyses 
(e.g. 0). 

Maximum-likelihood estimation is closely related to the optimal quadratic 
estimators (SUIIS- These take the form 

d * = \ E F m(i'R') tTSLCe ^ ssT ~ N)C- 1 aS/9C/?' C" 1 ], (33) 

I'R' 

where the Fisher matrix is FnfyQ'R') = tra.ce[C^ 1 dS/dC( i C~ l dS/dC^}/2. 
The quadratic estimator is unbiased for any (invertible) matrix C, but the 
variance is minimised if C is truly the covariance of the data. Of course, this 
is not known a priori so we instead construct an approximation to C from 
some reasonable guess for the power spectra. It can be shown that iterating 
the quadratic estimator (i.e. constructing C for the next estimate from the 
current power spectrum estimates Of 1 ) is equivalent to performing a Newton- 
Raphson maximisation of the likelihood, with the curvature approximated 
by the Fisher matrix |75| . The optimal quadratic estimator has the same 
computational problems as maximum likelihood methods. 

Heuristically-weighted, quadratic estimates 

The computational complexity of the optimal quadratic estimator ll-i-il) arises 
from the non-local weighting of the data by the inverse covariance matrix 
C _1 . Heuristically-weighted, quadratic methods ED IZ3 EH E2] circumvent 
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this problem by adopting a local weighting scheme. The use of fast-spherical 
transforms, as implemented in e.g. the HEALPIX package [53], reduce the 
operations count to O(N^) which is fast enough to be used in high- volume 
simulations. 

We shall concentrate here on the polarization auto-correlations only. In- 
cluding the cross-correlation between temperature and polarization is straight- 
forward |251 I82| . Assuming we are analysing Stokes maps, we adopt a real 
weighting w(n) that is zero where we have no data. For the case of temper- 
ature anisotropies, weighting with the inverse noise variance on scales where 
the noise dominates the signal, and uniformly where the signal dominates, can 
be shown to be close to optimal (TSj ■ Care should be taken in generalising this 
reasoning to low-amplitude B-modc polarization because of the additional 
complication of E-B mixing. We extract the electric and magnetic multipoles 
of the weighted polarization field with spherical transforms, 

E lm ±\B lm = J w{h){QTiU)(n) ±2 Y l * m (n)dn , 

= ^ ±2l(lm)(lm)'{E(i m y ±lB(i m y) , (34) 
(2m)' 

where the Hermitian coupling matrices are 



±il{lm){lm)> = J w(n) ± 2Y ilm y(n) ±2 Y l * n (n)dh , (35) 

and can be expressed in terms of the (scalar) multipoles, wi m , of the weight 
function, and the Wigner 3j symbols |82|. We have left the effect of smoothing 
by the instrument beam implicit in l|34|) ; the factors 2 Wi can be absorbed into 
the multipoles Ei m and B[ m , and, later, 2 Wj i into the power spectra Cf and 
Cf . We compress the pseudo- multipoles into pseudo-Ci s according to 



Cl =21 



^]T|^ m | 2 , Cf s-J—^lB^f . (36) 



These quantities are equivalent to the s T C 1 (dS/dCf)Q 1 s that appear in 
the optimal quadratic estimator (|33|l if we replace C _1 by heuristic weights. 



E-B mixing revisited 

We make a brief aside to note the relation between the pseudo-multipoles and 
the issue of separating E and B modes with incomplete sky coverage. The 
pseudo-multipoles mix E and B modes according to 

Elm = ^ { + I(lrn)(lrn)'E(i rn y + l-I(lrn)(lrn)> B(i m y) , (37) 
(lm)' 

Bim — ^ (+I(lm)(lm)'B(l m y — i-I(i m )(i m yE(i m y) , (38) 
(lm)' 
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where we have defined the Hermitian ±I(i m )(i m y t° be 

+ I{lm){lm)' = 2 (+2l(lm)(lm)' + -2-tym) (lm)' ) 
-I(lm)(lm)' = 2 (+2-tym)(im)' — -2^m(im)') = 

The mixing matrix -In m )n m y controls the mixing of E modes into B\ mi and 
vice versa. For the case of uniform weighting over the full sky, the mixing 
matrix vanishes and + In m \n m y cx 8u>8 mm >. For partial sky coverage, if w(n) 
is unity inside the observed region but zero elsewhere, integrating by parts 
in Q40(l shows that —In m sn m y reduces to a line integral around the bound- 
ary 



i-I {lm){lm y = ^ j-^jl <f (VbY^Y^f, - Y? m V h Y$$) dl a . (41) 

If we contract the Bi m with the multipoles Wj* of any scalar function 
W(n) — J2i m Wi m Yi m (n), such that W and \7 a W vanish on the bound- 
ary of the observed region, we see that J2im Wj* m -Iu m \(t m y = and we 
will have projected out the iS-mode contamination from the Bi m . Since the 
construction X^zm ^im^ lm 1S equivalent to integrating the contraction of 
Wg h * — v / 2e c (oVf ) )V c T^* and V a b over the observed region, isolating those 
functions in harmonic space that are orthogonal to the range of ~I(i m )(im)' is 
equivalent to finding elements of the tensor basis of the unambiguous E and 
B modes over the observed region (see Sect. 16.4^ . The harmonic-space method 
for E-B separation developed in [72] uses singular value techniques to project 
out the range of the mixing matrix -In m )u m y • 

For a more general weight function w(n), if we assume that it smoothly 
apodizes the edges of a connected region of linear dimension ~ i?, the relative 
sizes of -I and + I at scale Z _1 are in the ratio 1/(IR) |84|. It follows that 
the geometric effect of mode mixing is suppressed on scales small compared 
to the survey size. 

Mean values of the pseudo-Cjs 

Returning to the pseudo-C/s, their mean values contain two contributions, 
one from the covariance of the instrument noise (and foreground residuals) 
and one from the CMB power spectra. The former can be estimated by us- 
ing a large number of simulated observations: ~ 10 4 observations should give 
~ 1% accuracy. The mean level in the simulations can then be subtracted 
from the observed Ci to remove the noise bias in the mean. From now on 
we shall assume that the pseudo-C;s have been corrected in this way. Alter- 
natively, if only correlations between different instrument channels are used 
in the construction of the pseudo-C;s, the noise bias will vanish if the noise 



= / vY^^Y^dh, (39) 
J ^ r : )ah Y^dh .(40) 
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Fig. 8. Left: power spectra for B (top) and E (bottom; solid lines) compared to 
the mean pseudo-C;s (dotted lines) and the recovered power spectra (dashed lines). 
Right: the bottom panel shows representative matrices P ;i ' (solid lines) and M n i 
(dotted lines); the top panel shows the pseudo-inverses P^, 1 (solid lines) and M^} 
(dotted lines). Note that M u i and M^} have been multiplied by a factor of 100 
for clarity. The weight function applied to the map is uniform inside a circle of 
10° radius, with cosine apodization out to 15°. To obtain the pseudo-inverses, a 
Gaussian apodization of 4° HWHM is applied to the correlation functions. 



is uncorrelated between channels [S3]. The means of the pseudo-CjS are then 
linearly related to the true power spectra: 

(Cf ) = E P »' °v + Mw Cjf , = £ M W C# + Pw C» . (42) 

V V 

Here, the matrices are P w = J2 mm > \+I(lm)(lm)' | 2 /( 2 ^ + 1), with a similar 
expression for Mu>. They depend on the weight function w(n) only through 
the rotational invariant wi = J2 m \ w im\ 2 / (21 + 1), and can be expressed in 
terms of 3j symbols as 

{P/M) w = 2 I^1 £(2L + l)w L [l ± (-!)'+«'+*] ( l _ 2 ^)\ (43) 

For uniform weighting over the full sky, only wq is non-zero and so Mw van- 
ishes as it must (the 3j symbol forces I = I' and hence I + V + L to be even). 
For weight functions that smoothly apodize observations over a region of lin- 
ear size ~ R, the weight function is almost band- limited to l max ~ 1/R. In 
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this case the Pw and Mw matrices are close to band diagonal, and, for spec- 
tra that are smooth on the scale of Br 1 , the means of the pseudo-C;S are 
approximately 

(Cf) « Cf J2 Pw + Cf J2 M W , (Cf) « Cf J2 M W + Cf J2 Pw ■ 

V l' V V 

(44) 

The relative size of the mixing of E- and B-mode power in the pseudo-C;s is 
thus set by the ratio of the normalisations J2i> Mw to J2i> Pw ■ For I 3> 1/-R, 
it can be shown that [HI] 

E F «' w i /™ 2 («)dn, (45) 

E^'-il^/^Wdn, (46) 

so that mixing of power is suppressed by 1/(IR) 2 . Given that it is expected 
that Cf 3> Cf on all scales, {Cf) is essentially a scaled version of the true 
spectrum, but this need not be the case for Cf . These ideas are illustrated in 
Fig. |H1 which show typical rows of the matrices Py and Mu>, and the mean 
pseudo-Cjs that they give rise to, for observations covering a circle of radius 
15° with the last 5° apodized with a cosine function. The mean of Cf is seen 
to be a scaled version of the true power spectrum except on large scales, but 
this is not the case for Cf. 

'Unbiased' estimators 

The pseudo-C;s can be used directly for parameter estimation 76 , but it is 
often desirable to deconvolve the geometric effects of the weight function so 
we have unbiased estimates of the underlying power spectra. For example, for 
presentation purposes it is useful if we can plot the data in a form where they 
can be easily compared with true power spectra. 

For observations covering a significant fraction of the sky, deconvolution 
is straightforward. The matrices Pw ± Mw will then be invertible and we can 
form unbiased estimates of the power spectra according to 

Cf ± Cf = £(P ± M ) u }{Cf ± Cf) . (47) 
v 

An equivalent formulation of this approach can be made in terms of correlation 
functions |5S1|7S]. Consider estimating the correlation functions by aver- 
aging over pairs of pixels separated by /3, with pairs weighted by w(ni)w(n 2 )- 
It turns out that the data enters the estimators only through its pseudo-CjS. 
The estimators £±(/3) for ((Q =F iU)(hi)(Q + iJ7)(n 2 )), where hi ■ n 2 = cos/3, 
can be expressed as fZ5\ 
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2^ wiPi{cos (3) 

which also provides an efficient O(N^) way of estimating the correlation 
functions. If enough sky is available to estimate the correlation functions for 
all /3, then we can recover unbiased estimates of the power spectra by integral 
transforms. Using the orthogonality of the reduced Wigner functions, we have 

Cf ± Cf = 2tt J + £±03)4 ±2 03) dcos/3 . (49) 

The integration can be performed essentially exactly with Gauss-Legendre 
integration, and a significant feature of estimating the correlation functions 
via the pseudo-C;s is that they can be evaluated for any (3 up to the maximum 
pixel separation in the map. This avoids the need to re-sample the correlation 
functions onto the points required for the quadrature. This method has re- 
cently been applied to estimate the power spectra from the polarized channel 
of Archeops |66| . Aside from discretisation issues, this correlation function 
method is equivalent to performing the direct inversion in 14711 . 

For small surveys we will not have pixels separated by all angles < 180°, 
in which case we cannot estimate the correlation functions for all (3. The 
inversion in l)49|l is then not possible. This is equivalent to saying that the 
survey does not have the spectral resolution to perform the matrix inversions 
in l(4Tjl . Various schemes can be adopted to regularise the inversion, such as 
the use of bandpowers JZZJ- Alternatively, we can settle for pseudo-inverses 
Pj^ and M/7' 1 ' constructing estimators, e.g. 

CF = J2 P w^+M^, (50) 
v 

that have a simpler relation to the true power spectra than do the pseudo-Cjs. 
One desirable property we might enforce on the pseudo-inverses is that E-B 
mixing is removed in the mean, i.e. (Cf 3 ) has no contribution from Cf . A 
convenient way to find such a pseudo- inverse is given in |25| . which builds 
on earlier work in the context of cosmic shear analysis (SS] ■ We make use of 
the correlation functions estimates l|4*5|> in the angular range (0,/3 max ), where 
/?max is the maximum separation of pixels in the Stokes maps. It is shown 
in [25J th & t the estimators £(/3) ± £-(/?), which are linear in the pseudo-C;s, 
contain only E and £?-mode power, respectively, in the mean. The function 
£(/?) is obtained in the range (0, /3 ma x) by quadrature of in the same range: 

m = i+tf) + , 1 f £+(/?') sec 4 (072) dcos/3' 
sin {(3/2) J cosl3 

2(2 + cos/3) f 1 ta^/2) 

?+(P ) :— 5; — dcos£< , (51) 



sin 4 (/3/2) 7cos/3 ' sin/3' 
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and satisfies 

l(i W ±i_ m =j2^cf /B 4_ 2 (p) . (52) 
i 

We can recover estimates Cf and Cf that are linear combinations of the true 
power spectra, but with no mixing of E and B in the mean, by an apodized 
integral transform: 

Cf ,B = 2tt / ±[|(/3) ±t(/3)]4- 2 (/3)/(/3)dcos/3 , (53) 

J cos /3 ma * Z 

where /(/3) is chosen to apodize the discontinuity at /3 max - In this manner, 
we construct pseudo-inverses that deconvolve the pseudo-CzS up to a window 
function _ 2 Kw , i.e. {Cf') = J2 t , -iK\\iCft and similarly for Cf . The window 
function is determined solely by the apodizing function f{0) and is given by 

21' 4- 1 f 1 i 
-*K lv = —I-\ /(/?)4„ 2 (/3)4_ 2 (/3)dcos/3. (54) 

Oscillations in the window function can be suppressed with a careful choice of 
apodizing function, but the minimum achievable width is determined by the 
survey size l//3 max . An example of the construction of 'unbiased estimates' 
with this route is given in Fig. [S] 

Covariance of the pseudo-C;S 

The covariance matrix of the pseudo-C;s can be used to quantify the errors 
on any heuristically- weighted quadratic estimate of the power spectra. Having 
accurate error information is clearly essential to establish reliable constraints 
on cosmological models. For Gaussian fields, computing the pseudo-C; covari- 
ance directly is simple in principle, but the computational costs are prohibitive 
at high I. This is particularly problematic given that the sample covariance 
(i.e. that due to the random nature of the CMB fields) depends on the power 
spectra, and so should ideally be recomputed at each point in parameter space 
when constraining models. This model dependence also makes computing the 
sample variance 'blindly' from a large-suite of simulations unattractive. For- 
tunately, for the temperature anisotropies the sample covariance can be accu- 
rately approximated with analytic methods for I large compared to the inverse 
of the survey size |79j . Extending this calculation to polarization is difficult 
because the effect of E-B mixing can be the dominant source of variance in 
Cf on a wide range of scales for small surveys. However, useful analytic ap- 
proximations do now exist |84| . For an accurate calculation on all scales, one 
can combine the analytical calculations at high I (possibly calibrated from a 
small number of simulations if required) with a direct calculation at reduced 
resolution at low I [ZHj- (See also [HIj for an alternative prescription.) The noise 
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Fig. 9. Sample variance errors on the recovered Cf using the estimator <I53|I in the 
null hypothesis r = 0. The survey area and weight function are the same as in Fig- [HI 
The error boxes are the contribution purely from Cf (i.e. the lens-induced B modes) 
to the one-sigma errors on flat bandpower estimates with a Al — 70. These are thus 
representative of the errors that would be obtained with the pseudo-C; method if it 
were applied to maps from which the E modes had been removed (see Sect. 16.41 but 
here we have ignored the loss of information to the ambiguous modes). The error 
bars are the contribution purely from Cf and arise from E-B mixing. The dashed 
and solid lines are crude 'rule-of-thumb' approximations to the errors; see |84| for 
full details. 



contribution (and any residual errors from foreground removal) is generally 
best handled with simulations. 

An important issue for heuristically-weighted quadratic estimators is that 
the weighting scheme adopted generally makes no attempt to separate E 
and B modes coherently. While we can remove E-B mixing from the power 
spectra in the mean, in any realisation Cf need not vanish even if there are no 
-B-modes present. This cross-contribution to the sample variance of pseudo- 
C'i estimators makes them inaccurate for high-sensitivity, small-area B-mode 
surveys. This is illustrated in Fig. EI where it is shown that the contribution to 
the sample variance of Cf due to leakage of E modes is well in excess of that 
due to lens-induced B modes. (We have assumed there are no gravitational 
waves present, so lensing is the only source of B modes.) For the 15°-radius 
survey considered, with no instrument noise, this excess variance increases 
the error on the tensor-to-scalar ratio r in the null hypothesis to Ar = 0.15 
(3<r). One way to remove this problem is first to isolate the B modes with the 
methods described in Sect. 16.41 Ignoring the loss in this process, we find that 
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the limit on r is then reduced to 3.0 x 10 4 ; this arises solely from the sample 
variance of the lens-induced B modes. 

6.6 Non-Gaussianity and lensing reconstruction 

The effect of lensing on the CMB will probably be seen first in the angu- 
lar power spectra (e.g. |35l I88| and references therein), or, indirectly, via 
cross-correlation with large-scale structure |89| . As with gravitational- wave 
searches, the most promising place to look from the standpoint of cosmic 
variance is the i?-mode spectrum since, on small scales, lensing may be the 
dominant signal. Lensing is sensitive to the evolution of the gravitational po- 
tential at low redshift and can thus can be used to probe parameters, such as 
neutrino masses and the properties of dark energy, to which the primary CMB 
spectra are insensitive [HE]- At the sensitivity of the Planck experiment, the 
effect of lensing on the observed power spectra must be included to avoid bias 
when estimating parameters such as the baryon density. As thermal noise lev- 
els improve, and the sample variance of lens-induced -B-modes comes to dom- 
inate the error budget on Cf ' , corrections due to the non-Gaussianity from 
lensing will need to be included in the Cf errors and their covariance |74j . 

Looking further ahead, a number of authors have proposed reconstructing 
the lensing deflection field directly from the CMB by exploiting the non- 
Gaussianity of the lensing effect (see e.g. [5U1 EHI and references therein). 
This allows one to extract more information from the observed (lensed) CMB 
fields than is possible with an analysis of their power spectra alone, as well 
as providing a coherent reconstruction of the projected mass distribution. 
In essence, the reconstruction methods work by using the locally-anisotropic 
effect of lensing shear on the statistics of the CMB fluctuations within a 
region over which the deflections are coherent (~ 60°). The fact that we 
do not know the unlensed CMB fields, but only their statistical properties, 
introduces a noise in the reconstruction analogous to the effect of the scatter 
in intrinsic ellipticities of background galaxies in cosmic shear analyses. CMB 
polarization is particularly useful here |921 1931 I94j : it has more power on 
small scales than the temperature anisotropics and so allows reconstruction 
of the deflection to smaller scales. It was shown in that local correlations 
between E and the lcns-induccd B modes, estimated with simple quadratic 
estimators, reconstruct the deflection field with the highest signal-to-noise on 
all scales once thermal noise permits imaging of the lens-induced B modes 
on small scales. This requires both high sensitivity (~ 1 fiK in an arcmin 
pixel), and also high angular resolution (a few arcmin). The E-B estimator 
is particularly powerful since there is no contribution to its sample variance 
from chance correlations in the unlensed fields since there are no unlensed 
small-scale B modes. 

Recently, it has been argued that the relatively larger non-Gaussianity of 
the small-scale (lensed) 5-modes, compared to the temperature anisotropies 
and E modes, allows further improvement in the lensing reconstruction if more 
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optimal techniques are employed |94| . Indeed, a simple counting argument 
suggests that, in the absence of thermal noise, the assumption that lensing 
is the only source of small-scale B modes allows perfect reconstruction of 
the deflection field (except possibly for a few degenerate configurations). In 
practice, errors in the theoretical modelling of lensing ultimately limit the 
reconstruction, but this effect is unimportant until thermal noise levels fall 
well below the very optimistic level of 0.25 /xK-arcmin. 

An important application of lensing reconstruction is to 'clean out' the 
large-scale B modes where any gravitational wave background is expected 
to contribute. If the lens-induced modes are treated simply as an additional 
source of Gaussian confusion, the limit on the detectability imposed by lensing 
(i.e. with no thermal noise or errors in foreground removal) is r > 4 x I CP 5 
at 3(7 for a full-sky survey |44| . This can be improved by more than a factor 
of 40 with optimal reconstruction methods corresponding to an energy 
scale of inflation V 1 ^ « I x 10 15 GeV. 



7 Conclusion 

Detections of CMB polarization are still at an early stage, but already we 
are beginning to see the promise of polarization data for constraining cosmo- 
logical models being realised. The electric-polarization power spectrum mea- 
surements from DASI and CBI jS] reveal acoustic oscillations with an 
amplitude and phase perfectly consistent with the best-fit adiabatic models 
to the temperature anisotropies. This is a non-trivial test on the dynamics of 
the photons and baryons around the time of recombination. Furthermore, the 
large-angle measurements of the temperature-polarization cross-correlation 
from WMAP [S] indicate a large optical depth to reionization and hence a 
complex ionization history. We can expect further rapid progress observa- 
tionally, with more accurate measurements of the power spectra of E-mode 
polarization, and its correlation with the temperature anisotropies, expected 
shortly from a number of ground and balloon-borne experiments. These exper- 
iments will also greatly increase our knowledge of polarized astrophysical fore- 
grounds at CMB frequencies. Instruments are currently being commissioned 
that should have the sensitivity to detect the power of -B-mode polarization 
induced by weak gravitational lensing |95j . and, already, several groups are 
working on a new generation of polarimeters with the ambition of detecting 
gravitational-waves and reconstructing the projected mass distribution from 
CMB polarization observations. 

The success of these programmes will depend critically on many complex 
data analysis steps. We have attempted to summarise here some of the generic 
parts of a polarization data analysis pipeline, but, inevitably, have had to leave 
out many topics that are more instrument-specific. Important omissions in- 
clude calibration of the instrument, cleaning and other low-level reductions of 
time-stream data, noise estimation, propagation of errors, and the broad topic 
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of statistical accounting for non-ideal instrument effects. Given the exquisite 
control of systematic effects - both instrumental and astrophysical - that 
searching for sub-/j,K signals demands these omissions will almost cer- 
tainly prove to be the most critical steps. 
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